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ABSTRACT 

2M1938+4603 (kic 9472174) displays a spectacular light curve dominated by a strong 
reflection effect and rather shallow, grazing eclipses. The orbital period is 0.126 days, 
the second longest period yet found for an eclipsing sdB+dM, but still close to the 
] minimum 0.1-d period among such systems. The phase- folded light curve was used 

Ti^lj- ■ to detrend the orbital effects from the dataset, and the resulting amplitude spectrum 

shows a rich collection of pulsation peaks spanning frequencies from ~50 to 4500 |xHz. 
The presence of a complex pulsation spectrum in both the p-mode and the g-mode 
regions has never been seen before in a compact pulsator. 

Eclipsing sdB+dM stars are very rare, with only seven systems known and only 
one with a pulsating primary. Pulsating stars in eclipsing binaries are especially im- 
portant since they permit masses derived from seismological model fits to be cross 
checked with orbital mass constraints. We present a first analysis of this star based 
■ on the Kepler 9.7-day commissioning light curve and extensive ground-based pho- 

' tometry and spectroscopy that allow us to set useful bounds on the system param- 

eters. We derive a radial- velocity amplitude if i = 65.7 ±0.6 km s^^, inclination angle 
i = 69.45 ±0.20 °, and find that the masses of the components are Mi = 0.48 ±0.03 
and Ma = 0.12 ±0.01 Mq. 

Key words: subdwarfs - binaries: close - binaries: eclipsing - stars: variables: general 
- stars: individual: 2M1938+4603 
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1 INTRODUCTION 

The subdwarf B (sdB) stars are known to be core helium 
burning stars with extremely t hin (Menv ^ -01 Mp) inert 
hydrogen dominated envelopes ([Hebeif 1 1986'). This places 
them on an extension to the classical horizontal branch. 
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known as the extreme horizontal branch (EHB). In order 
to reach this configuration, almost the entire envelope must 
be stripped of! close to the tip of the red giant branch. There 
are several binary scenarios that are capable of accomplish- 
ing this, such as common-envelope ejection (CEE), stable 
Roche-l obe ove rflow and merger of two helium-core white 
dwarfs (|Han et al. 2002). A small minority of sdB stars in 
the field are found to have close M-dwarf companions (see 



© 2010 RAS 



L2 R. H. 0stensen et al. 



iFor et alJlioiol . for recent results and a review), and seven 
such systems are eclipsing. The eclipsing systems have been 
monitored over long time-bases in order to detect low-mass 
companions from precise m easurements of the eclipse tim- 
ings with the 0-C method (|Kilkennv et alJboOSl : iLee et alj 
|2009| ). 



One of the eclipsing sdB-l-dM systems, NYVir, 
has a pulsat i ng p rimary of the V361 Hya class 
jKilkennv et al.1 Il998l'). The V 361Hya stars were dis- 



covered bv iKilkennv et al] (|l997l ) and are characterised by 
rapid pulsations, typically in the period range between two 
and five minutes. They are known to be pressure (p-)mode 
pulsators excited by the k mechanism, driven primarily 
by an iron opacity bump in the envelope (|Charpinet et al.] 
I997I). A se cond class of sdB pulsators was reported by 
Green et al ] ([2003). These stars, now known as V1093Her 
stars, show pulsations with much longer periods (~1 h) than 
the V361 Hya stars, and their temperatures are lower. These 
pulsations can be described in term s of gravity ((7- ) mode s 
excited by the same k mechanism (|Fontauie et alj boOSh . 
With the discovery of l ong-period pulsati ons in a known 
rapid pulsator, DWLvn. lSchuh et al.l ()2006l ) established the 
existence of hybrid sdB pulsators. A broad review of ho t 
subdwarf stars in general can be found in iHebeJ l|2009l l. 
and a review of aster oseismology and evolution of the EHB 
stars can be found in l0stenserj (|2009|). For the most recent 
pulsator discoveries see l0stensen et al] l|2010al ). 
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Figure 1. The first half day of the Kepler light curve. The upper 
panel shows the raw data, which is dominated by the strong re- 
flection efl'ect, with both primary and secondary eclipses clearly 
visible. The middle panel shows the same chunk after detrending 
with the light curve folded on P, and the bottom panel shows the 
same after prewhitening the 55 most significant frequencies. 



The Kepler spacecraft was launched in March 2009, 
aiming to find Earth-sized planets fro m photometric obser - 
vations of a 105 square degree field (|Borucki et al.|[201oh . 
Kepler is also ideally suited for asteroseismological studies, 
and the Kepler Ast eroseismic Science Consortium (kasc, 
iGimiand et al 2010) manages this important aspect of the 
mission. The methods with which the compact pulsator can- 
didates were selected, together with analysis of the first half 
of the survey phase, are presented in l0stensen et aL I (|2010bl . 
Paper l) . The first resu lts on a V361 Hya star in the Kepler 
field are presented by iKawaler et al. I i|2010al ). the first re- 
sults on V1093H er and DWLyn pulsators are presented by 
iReed et all (|2010l ). and results on two VI 093 Her pulsators in 
sdB+dM reflection binaries are discussed in iKawaler et al.l 
l|2010bl ). Further analysis of Keple r data on sdB stars are 
found in IVan Grootel et ahl (|2010l ). who present the first 
detailed astero s eismic analysis of a V1093 Her star, and in 
iBloemen et al.l (|2010l ). who present a detailed analysis of 
the extraordinary binary light curve of the sdB-|-WD star 
KPD 1946-F4340. 



In this letter we present 2M1938-I-4603 (kic 9472174), 
an object first identified as an sdB star in the Kepler field 
durin g a survey of blue targ ets selected from 2MASS photom- 
etry ISkrutskie et al.ll2006l ). Follow-up photometry showed 
the presence of a reflection effect with shallow primary and 
secondary eclipses (Fig. [TJ . The amplitude of the reflection 
effect is very comparable to those observed in HWVir and 
NYVir, but the eclipses are much more shallow. After de- 
trending the Kepler light curve, we clearly detect low-level 
pulsations. 



2 OBSERVATIONS 

The discovery of the strong reflection effect with grazing 
eclipses was made by two of us (JBL & MM) during a pho- 
tometric run in June 2008. Here we present only the eclipse 
timings from the ground-based photometry, as the pulsa- 
tions are too complex and have too low amplitudes to be 
significant in those light curves. The times of thirteen pri- 
mary eclipses collected between June 2008 and May 2010 
(Table IXTj) give the following ephemeris: 

To = 2454640.864162 ± 0.000058 d 
P = 0.125765300 ± 0.000000021 d. 

We downloaded the Kepler QO light curve from the 
KASC archive, corrected the raw time-stamps to Barycentric 
Julian Date (BJD) according to the instructions (|Van Clevd 
[2009), and made minor corrections to the raw fiuxes to take 
out trends on time-scales longer than a day. The 77 consecu- 
tive eclipse times measured from the Kepler photometry are 
listed in Table [X2l and match the ground-based ephemeris 
well within the errors. 

In order to produce a useful Fourier Transform 
(FT) that shows the spectrum of low-level pulsations in 
2M1938-I-4603 among the extremely dominant orbital effects 
(Fig. [T] top panel) we first attempted to clean out a model 
light curve of the system, but we abandoned this approach 
since the FT of the residuals contain significant peaks at 
every orbital harmonic. The model light curve is unable to 
produce a satisfactory fit to the complicated irradiation ef- 
fect, at the exceptional precision of the Kepler photometry. 
Our model light curve does not account for radiative transfer 
through the heated face of the M-dwarf, which may account 
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Figure 2. The FT of the data after detrending the orbital and long term variations (upper part of each panel). The panels are 
labeled A, B, C, D, E and these also define the group designation used in Table IdT] The magenta arrows indicate the orbital frequency of 
farb = 92 |i.Hz and all its harmonics, which have been completely removed by subtracting the data folded forh- The green arrows indicate 
the Kepler long-cadence cycle, /ic and its harmonics, at which possible artefacts are known to be found. The lower part of each panel 
shows the FT of the data after prewhitening 55 frequencies, as marked with blue arrows. The blue line indicates four times the mean 
noise level at 24 |j.ma. 



for some of the discrepancies between the model and data. 
Until such issues are resolved the precise parameters of our 
model are subject to systematic uncertainties that could well 
be in excess of the statistical errors. Even if we managed to 
model all orbital effects in the light curve to the required 
precision, any pulsation peaks found in the residuals on an 
orbital harmonic frequency would still be suspicious. So in- 
stead we proceeded by folding the Kepler light curve on P, 
and then using the result to clean out all orbital effects from 
the light curve. Note that any stable pulsations that coin- 
cide with these periods within the resolution are also very 
effectively removed. 



2.1 Frequency determination 

After removing the orbital effects from the light curve, a 
spectacularly rich pulsation spectrum is revealed (upper half 
of the five panels in Fig. We have identified 55 frequen- 
cies between ^^50 and 4500 viHz that have amplitudes well 



above four times the mean level of 6 vimsQ computed from 
the whole FT between 500 yiRz and the Nyquist frequency, 
/Nyq = 8496 [iHz. These frequencies are listed in Table [Dl] 
Several more frequencies remain above the 24 jima limit, af- 
ter prewhitening the light curve with these 55 frequencies 
(lower half of the panels in Fig. [2)), especially in the low 
frequency region below 500 ji-Hz. However, we chose to con- 
strain ourselves to the most clearly resolved peaks in this 
first analysis, as much more Kepler data is underway, and 
will provide a tenfold increase in resolution and a threefold 
drop in the noise level after only the first three-month cy- 
cle of observations. In total we hope to obtain close to five 
years of almost continuous space based observations on this 
object, so we will limit ourselves to discussing only the most 
obvious features here. 

The FT shows several distinct groups of frequencies, 
and we number the frequencies in each group separately 

^ One micro-modulation amplitude (|^ma) is equivalent to a semi- 
amplitude of one part per million of the mean light level. 
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as A,B,C, D,E, as indicated in Fig. [2l The frequencies in 
group A are typical for the V1093Her stars, and we iden- 
tify and remove 20 peaks, with the lowest at 6 cr. There 
are clearly many more frequencies than this, in particular 
between 100 and 200 M-Hz, but the resolution is not good 
enough to justify pushing the limit further. There are no sig- 
nificant frequencies between 500 and 1000 \xRz, except one 
at 570 viHz which is just above the detection limit. Note also 
that this frequency is close to the long cadence readout cy- 
cle at /ic = 566.4 (j.Hz. In group B we find only six rather 
low-level peaks. Group C contains the strongest peak in the 
FT at fci = 2265.8 ^Hz with an amplitude of 423 M-rna. This 
frequency is almost exactly at 4/ic = 2265.4 viHz, but since 
we do not observe such strong artefacts in any other stars 
we are inclined to believe that this is a real pulsation mode. 
In all other short cadence Kepler datasets we have worked 
with, the strongest artefact peak is found at 8 or 9/ic. In 
this dataset 9/ic is not present at all, and 8/ic = 4531 (iHz 
is barely above the detection limit, while fci is 70-cr! The 
second strongest peak in the FT is also found right next to 
it, /c3 = 2200.1 [iHz with an amplitude of 233 jxma. Another 
rich set of peaks is found around 2700 H-Hz, so that group C 
contains 18 peaks altogether. Group d contains 9 clear peaks 
between 3500 and 3800 fi.Hz, and group E contains only two 
significant peaks located at 2 fa and 2/c4 = 8/ic. The ra- 
tio of some of the high-amplitude peaks in the short-period 
range, such as fm/fc4, ~ 1.638 and fog/fcs = 1.717, which 
are both close to the theoretical ratio of the radial funda- 
mental and the second overtone (1.667), and a stellar model 
for the Teff and logg of the primary does display pulsation 
periods for radial modes close to these frequencies. 

The frequency analysis and er ror estimations on th e 
parameters are done as outlined in iDegroote etlo] (|2009l ). 
Because some of the frequencies were closer together than 
the Rayleigh T^^ limit, and thus unresolved, we have per- 
formed additional analyses, where after each prewhitening 
stage, some of the frequencies, amplitudes and phases were 
refitted with a Levenberg-Marquardt nonlinear fitting rou- 
tine. To investigate the stability of the fit and the differences 
from the linear analysis, we first only updated the ampli- 
tudes, phases and frequencies belonging the last two identi- 
fied peaks, then updating the parameters belonging to the 
last 20 peaks, and finally using all identified frequencies. The 
differences between the amplitudes and frequencies were all 
within the derived error estimates, except for the frequencies 
around 3712 and 2593 ^.Hz. The nonlinear fitting algorithm 
did not produce consistent results at these frequencies, but 
future observations will allow us to resolve these structures. 
Below 200 (i-Hz, a dense forest of barely significant peaks 
is visible. Although the different methods give consistent 
results, more Kepler photometry will be necessary to fully 
resolve them. 
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Figure 3. Mass— radius diagram for 2M1938-I-4603 showing the 
regions permitted by the orbit solution (green curves) and by 
the spectroscopic gravity (red curves). The dotted, green curves 
show the error ranges that corresponds to a 3-(t error in K\ for 
the orbit solution. The dotted, red curves corresponds to the one 
(T formal error on log g for the gravity track. The green curves are 
labeled with the value for q = M\ / M2 that corresponds to each 
point on the curve, with + marking ticks of 0.01, increasing to 
the left. Note that changing K\ by Sct shifts the curve very little 
in the M /R plane, but the g- value changes considerably. Vertical 
lines show the mass range typical for canonical EHB stars formed 
through the CE channel. 

the exposures. The resulting spectra have a resolution of 
7? ~ 2150 over the range 3675 -4520 A (Bok), or i?~4200, 
4000 -4950 A (MMT). Radial velocities were derived by 
cross-correlating the individual continuum-removed spectra 
against a super-template using the iraf task FXCOR. The 
super-template for the lower resolution spectra was obtained 
by shifting the 66 individual spectra to the same velocity 
prior to median filtering into a single spectrum; the veloc- 
ity zero point of the resulting template was thus undeter- 
mined. The cross-correlation template for the 6 MMT spec- 
tra was constructed from 19 spectra previously obtained 
with the identical spectroscopic setup for other hot subd- 
warfs of known radial velocities, whose spectral abundance 
patterns closely match that of 2M1938-I-4603. 

A single sinusoidal fit was performed to all 72 veloci- 
ties as a function of orbital phase, weighted by the velocity 
errors and also including an additional term for the zero- 
point offset of the Bok velocities relative to the MMT ve- 
locities (Fig. IB1|) . The orbital period was fixed at the value 
of 0.12576530 d derived from the eclipse timings. The de- 
rived RV semi-amplitude is /(Ti = 65.7 ±0.6km s~^, with a 
systemic velocity 7 = 20.1 ±0.3 kms~^. 



2.2 Spectroscopy 

Spectroscopic observations to determine radial velocities for 
2M1938+4603 were undertaken primarily with the B&C 
spectrograph on the 2.3-m Bok Telescope on Kitt Peak. A 
few spectra were also obtained using the MMT Blue spec- 
trograph. In both cases, an 832/mm grating was employed 
in second order, and particular care was taken to main- 
tain precise centering of the star on the slit throughout 



3 SYSTEM PARAMETERS 

The effective temperature, surface gravity, and helium abun- 
dance of 2M1938-I-4603 were determined in the context of 
Paper I to be T^a = 29 564 ±106 K, log g = 5.425 ± 0.009 dex, 
and log (A'ho/A^h) = -2.36 ±0.06 dex, using metal blanketed 
LTE models. This temperature is close to the boundary re- 
gion between the p- and g-moAe pulsators, where hybrid 
DWLyn type pulsators have been found. 
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We modelled the light curve using grid elements cov- 
ering each star, accounting for tidal (ellipsoidal) distortion 
of each star and gravity and limb-darkening of the sdB as 
m iBloemen et~al] |2010l ). Irradiation of the M-dwarf was ac- 
counted for by summing the M-dwarf flux with the incident 
flux from the sdB point by point. Surface brightnesses were 
computed assuming black-body spectra at a single wave- 
length of 600 nm. The best fitting model is found to have the 
following parameters: inclination angle, i = 69.45(2)°, rela- 
tive radii, ri = 0.250(1), r2 = 0.177(1) in units of the orbital 
separation. The errors in parentheses are formal fitting er- 
rors that may underestimate the true errors perhaps as much 
as a factor 10, considering the discrepancies in the light curve 
fitting. The high precision on these parameters allows us to 
use the mass-radius relationships as derived from the or- 
bital parameters and from the surface gravity to constrain 
the mass and radius of the primary, as shown in Fig. [S] We 
clearly see that the permitted M /R given by the photomet- 
ric orbital parameters (P, i, ri and r2) and the spectroscopic 
Ki crosses the M /R given by the spectroscopic surface grav- 
ity at precisely the expected primary mass for a post-CE sdB 
star. The adjacent dotted lines indicate the errors (on the 
RV for the g-track, and on log g for the p-track since the er- 
rors on P, i and the relative radii are too small to matter). A 
mass for the primary of M\ = 0.48 ±0.03 M© can be deduced 
from the diagram. Using K\ and P we know the mass func- 
tion /(M)= 0.003695. With 0.48 Mq for the primary, and 
solving for the secondary, we get M2 =0.12 ±0.01 M©. 



4 DISCUSSION AND OUTLOOK 

We have presented the discovery of a new eclipsing sdB with 
an M-dwarf companion, and demonstrated that the primary 
is an unusual pulsator with frequencies spanning the range 
from 50 to 4400 (j.Hz. The unprecedented precision of the Ke- 
pler observations makes a direct comparison to other known 
eclipsing sdB+dM systems rather difficult, since such low- 
level pulsations as we find in 2M1938+4603 would be un- 
detectable from the ground. Nevertheless, it is an amazing 
stroke of luck to find an eclipsing sdB+dM system within the 
narrow confines of the Kepler field, especially one as bright 
as 12th magnitude, considering the fact that only seven such 
systems are known among the thousands of subdwarfs that 
have been surveyed for variability to date. 

In future papers we will produce a more realistic model 
that can handle the irradiation effect of the M-dwarf to im- 
prove the binary light curve model, and produce asteroseis- 
mic models that can reproduce the frequencies observed in 
this star. When longer time-series of Kepler photometry be- 
come available, we will be able to fully resolve the densely 
packed amplitude spectrum that we see in the g-mode re- 
gion between 100 and 200 [iHz. This will enable us to pro- 
duce an asteroseismic model reproducing all the frequencies 
in this complex rotatin g pulsator, as was done for NY Vir by 
ICharpinet et al.l (|2008^. Since 2M1938+4603 has p-modes in 
addition to p-modes we should also be able to constrain the 
deeper structure of the sdB core, and hopefully also establish 
the mass of its progenitor, before the envelope was ejected 
by the companion. With five years of precise eclipse timings 
we may detect planetary companions from variations in the 
O-C diagram, if they orbit this compact binary with peri- 



ods shorter than the Kepler Mission, and the low-amplitude 
pulsations are found to be stable on such time-scales. 
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Table Al. Times of primary minima from ground based obser- 
vations, and O-C values for To = 54640.864162, P = 0.125765300. 



Epoch BJD CTBjD O-C 



Q 


54640, 


.86420 


0, 




0, 


00004 


595 


54715, 


.69440 


0, 


.00004 


-0, 


.00012 


714 


54730, 
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0, 


.00004 


0, 


.00001 
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.93960 


0, 
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Table A2. Times of primary minima from Kepler observations, 
and O-C values for To = 54640.864162, P = 0.125765300. 



Epoch 


BJD 


0"BJD 


O-C 
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0, 
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-0, 


.00000 
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.06262 


0, 


.00001 


0, 


.00001 
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Table A3. Times of primary minima from Kepler observations, 
cont. 



Epoch BJD CTBJD O-C 



2547 


54961, 


,18842 


0, 
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,31416 


0, 


.00001 


0, 


.00001 


2549 


54961, 


.43995 


0, 


.00001 


0, 


.00004 


2550 


54961, 


.56569 


0, 


.00001 


0, 


.00001 


2551 


54961, 


.69145 


0, 


.00001 


0, 


.00001 


2552 


54961, 


.81719 


0, 


.00001 


-0, 


.00002 


2553 


54961, 


.94298 


0, 


.00001 


0, 


.00001 


2554 


54962, 


.06873 


0, 


.00001 


-0, 


.00001 


2555 


54962, 


.19452 


0, 


.00001 


0, 


.00002 


2556 


54962, 


.32026 


0, 


.00001 


-0, 


.00001 


2557 


54962, 


.44603 


0, 


.00001 


-0, 


.00000 


2558 


54962, 


.57180 


0, 


.00001 


0, 


.00000 


2559 


54962, 


.69755 


0, 


.00001 


-0, 


.00001 


2560 


54962, 


.82334 


0, 


.00001 


0, 


.00001 


2561 


54962, 


.94909 


0, 


.00001 


-0, 


.00001 


2562 


54963, 


.07486 


0, 


.00001 


-0, 


.00000 


2563 


54963, 


.20060 


0, 


.00001 
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orbital phase 

Figure Bl. The spectroscopic RVs. Open points are from the 
Bok Telescope, filled dots from the MMT. 
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Orbital phase 

Figure CI. Best fit model light curve and the folded Kepler 
light curve (top), and residuals (bottom). 
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2M1938+4603: An eclipsing sdBV+dM observed with Kepler L9 



APPENDIX D: FREQUENCY LIST 

The two pairs of frequencies marked with ^ are unresolved. 
Prewhitening one frequency with the measured amplitude 
leaves a significant residual peak, but fitting two frequen- 
cies simultaneously produces possibly exaggerated ampli- 
tudes and high errors on the frequencies. The amplitude 
given is therefore only stated as the observed amplitude of 
the combined peak. 



[h!] 

Table Dl. Clearly resolved frequencies, periods, and amplitudes 
for 2M1938-I-4603. 
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